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Abstract 

We study the discovery potential of the CMS detector for the scalar leptoquark pair production 
at the LHC. Present and future exclusion limits are considered. We find that the maximal leptoquark 
mass reach is m ~ 1.47 TeV for the branching ratio of Br(LQ; — ► 1 q) = 100 %, while for Br(LQ; — > 1 
q) = 50 % the upper limit is 1.2 TeV for an integrated luminosity of 100 fb~ . We obtain comparable 
results for electron and muon-type leptoquarks. The pileup effect at high luminosity is discussed. 
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1 Introduction and phenomenology 



Leptoquark particles appear in extensions of the Standard Model like Grand Unified Theories, composite 
models and R-parity violating Supersymmetry. They carry both lepton and baryon number. These 
particles mediate transition between quark and lepton . Leptoquarks can be scalar (Spin= 0) or vector 
(Spin=l) and couple to leptons via a Yukawa type coupling. These particles conserve lepton and baryon 
number in order to avoid a fast proton decay. One can classify the different leptoquark through their 
spin, weak isospin and fermion number, see details in [Q, Q]. 

In the following we consider only pair production of scalar leptoquarks: almost no dependence on the 
A coupling is expected (coupling to lepton and quark). The leptoquark resonance width can be defined 
as (El): 



T ~ 350MeV 
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We obtain r ~ 2 GeV for A = y 'Ana and Mlq — 1200 GeV. As we will see later the theoretical width 
is overwhelmed by the calorimeter resolution. If we assume that leptoquarks have either left or right 
couplings but not both, the branching fraction is 0, 0.5 or 1 for symmetry reason. We consider both 
electron and muon leptoquark type with Br(LQ e — ► e q) = 1 or Br(LQ e — > c q) = 0.5 and Br(LQ M — * fx 
q) = 1 or Br(LQ„ -> n q) = 0.5. 

Fig.l shows the leading-order diagrams for the leptoquark pair production at LHC. One can notice 
that except for the diagram of the type qq — > LQLQ via e-exchange, all the other processes do not depend 
on the A coupling (only vertex g — LQ — LQ). This is very important as the cross section then depend 
mainly on the mass of the leptoquark. 
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Figure 1: Leptoquark pair production diagrams 



We take the total cross section for the leptoquark pair production from Q (Next-to-Leading Order 
calculations (NLO) using CTEQ4M structure function). The difference between LO (as in PYTHIA {§) 
and NLO can reach 50% depending on the leptoquark mass. Thus NLO cross section is taken for our 
analysis. 



2 Exclusion limit 

The current limit on the leptoquark originates from HERA and Tevatron ||, |(| analysis. At the Tevatron 
the exclusion limit for NLO cross section is investigated by both CDF and DO jjj. Combining both 
results one obtains : Mf Q lar > 246 GeV. We can expect at y/s — 2 TeV with L = 1 fb _1 for the Tevatron 
(D0+CDF) a limit of Mf c £ lar > 350 GeV. 

At LEP the exclusion limit was investigated, e.g., by OPAL using the fact that the leptoquark enhances 
the cross section of e + e~ — > qq in t/u channel exchange ||. A preliminary study carried out in CMS 



|H showed that the mass reach is about 1.5 TeV 
LHC performance of leptoquark detection: O, n3, 



Several studies have been done to evaluate the 
The study presented in this note includes a more 



accurate detector performance description (as will be explained in next section), takes into account all 
sources of background and includes superposition of events at high luminosity (so called pileup) . 
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3 Simulation and CMSJET/CMSIM comparison 

Fig. 2 shows the full GEANT simulation of a 900 GeV leptoquark pair decay into electron (positron) and 
u (u) quark jet (pp — » LQ e LQ e — > e~e + uu) in the CMS detector. 




Figure 2: Full GEANT simulation of pp -> LQ e LQ e -> e-e+uu for m LQo = 900 GeV. 

The GMSJET fast Monte-Carlo package [[t4| is used to model the CMS detector response due to 
the huge sample of events to be processed (about 20 millions at generator level and half a million at 
the simulation level). The fast simulation incorporates the full granularity of the calorimeter with an 
energy smearing according to GEANT simulation and parameterization of the longitudinal and lateral 
profiles of showers with some cracks description (transition regions of calorimeter: endcap-barrel, endcap- 
forward calorimeter), charged particles and the muon resolution is parameterized as function of px and 
r\. Here we illustrate the performance of the fast simulation with some distributions in Fig. 3, where the 
1 TeV jet response is produced by both fast simulation with CMSJET and full detector simulation with 
CMSIM We compare the energy resolution, transverse energy distribution (i. e. the ratio of the 

energy contained in a cone of R = y/rf + <j) 2 value to the total energy) and the fraction of the jet energy 
deposited in electromagnetic calorimeter. 

4 Background 

The background processes considered in the following are required to have 2 isolated high-py leptons and 
at least 2 energetic jets. Due to the large QCD cross section, semileptonic decays of heavy flavour are 
also considered requiring isolated high-p^ leptons (coming from jets). 

In Tab.l we give the cross section for some processes where the production of 2 leptons and 2 jets is 
forced in PYTHIA. For some processes with a very large total cross section the generated samples are 
subdivided in several pr intervals. 
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Figure 3: Full simulation versus fast simulation for 1 TeV jets 



Due to the large QCD cross section a special treatment is applied for those events. After parton 
generation and hadronization each meson with sufficient value of pr is decayed 1000 times in order to 
select events where at least 2 high-px leptons (pr > 40 GeV) of same flavour and opposite sign (SF OS) 
are produced (see Tab. 2). Once passes this requirement the QCD event is processed through CMSJET 
simulation. This procedure increases the generated statistics by a factor 1000. After this simulation no 
QCD event survives the first level selection: 2 isolated high-p^ leptons (pt > 40 GeV in |?7| < 2.4) of 
opposite sign and same flavour and at least 2 jets ( E<r > 40 GeV in |?7| < 4.5). 



Process 


ZZ(lljj) 


ZW(11 jj) 


ti(l v j 


Zjet 
(jot > 20 GeV) 


Wbt(l v 1 v jj) 
(pr > 20 GeV) 


WW(1 v 1 v) 


a (fb) 


600 


596 


7 • 10 a 


1.25 • 10 5 


10 3 


827 


Number of events 


6- 10 4 


6 • 10 4 


7 • 10 6 


1.25- 10 7 


10 5 


8.27- 10 4 


Nbr of gen. evts 


1.8 • 10 5 


1.8 • 10 5 


10 6 


8.4- 10 B 


2 • 10 5 


10 5 



Table 1: Background cross section and number of events for L = 10° pb 



5 Selection variables 

The selection is based on several basic variables with different cut values to optimize the signal over 
background ratio for various leptoquark masses. We require exactly 2 isolated leptons SF OS with pr > 
40 GeV and at least 2 jets with Et > 40 GeV (so called first level cuts). Fig. 4 shows some distributions 
after applying first level cuts for signal (leptoquark pair production of t^lq c — 900 GeV) and background. 
In order to reduce the background originating from Z decay in 2 leptons we apply the lepton pair mass 
cut: m;; > 100-200 GeV. The ti background is reduced by constraining the missing transverse energy to 
a low value : E™ ISS < 80-150 GeV. This reduces the amount of events where an energetic neutrino is 
produced. We also use the variable £ Ey which is the sum of the event ( cell by cell) transverse energy in 
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Figure 4: Distributions over selection variables (signal in light grey and background in dark grey, tulq^ 
= 900 GeV). 
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Kinematics 


a (fb) 


gen evts 


ratio 


100 GeV < PT < 200 GeV 


1.3 • 10 9 


1.1 • 10 9 


0.01 


200 GeV < pT < 400 GeV 


7.14- 10 7 


1.6 • 10 8 


0.02 


PT > 400 GeV 


2.8 • 10 6 


1.24- 10 8 


0.44 



Table 2: QCD cross section and number of events for L = 100 fb 1 



the calorimeter system down to |7j| — 4.5. The requirement on this value takes into account the pileup (see 
next section for details) which produces a shift of about 800 GeV, so we impose: SEr> 1500-2000 GeV. 
We also use the ratio of the last two variables to cope with the fact that the transverse missing energy 
tail in the case of balanced events ( no energetic neutrino ) originates mainly from a mismeasurement of 
high energy jets. Thus the variable E™ ss /E should be suitable to solve this problem as illustrated 
in Fig.4 . We require : E™ ss /E E T < 0.03-0.05. 

One of the most important requirement is related to the characteristics of pair production, therefore 
both pairs of the lepton-jet invariant mass has to be within a window: A my < 150-320 GeV. Finally, in 
order to obtain the best significance we search for a signal peak as an excess of events over exponentially 
falling background distribution with an optimised window: Am < 100-220 GeV. The procedure to obtain 
the best significance for several leptoquark mass values makes use of a loop over the previous cuts. Here we 
consider all the possible combinations between the 2 leptons and the jets (at least 2), resulting sometimes 
in more than 2 accepted lepton-jet pairs in one event. The same treatment is applied to the signal and 
to the background. 

A typical selection for m L Q = 1400 GeV is based on the following cuts: 

• 2 isolated leptons (SF OS) with p T > 40 GeV, \r]\ < 2.4; m u > 150 GeV 

• At least 2 jets with E T > 60 GeV , \r]\ < 4.5 

• Erp iss < 185 GeV, E E T > 1700 GeV, E^ ss /E E T < 0.04 

• A m tj < 310 GeV, the peak window Am = 210 GeV 

6 Effect of pileup on the observables 

The pileup originates from the fact that several protons can interact during the same bunch crossing. 
Assuming a total proton-proton cross section of 100 mb at the LHC (from a Regge-based extrapolation 
of low energy data) and taking into account the time interval between bunch crossing of 25 ns, we obtain 
an average of 25 interactions per bunch crossings (Poisson distribution) for an instantaneous luminosity 
of 10 34 cm^V 1 . In our simulation we superimpose these events produced with MSEL=2 by PYTHIA 
(it includes elastic scattering, single and double diffraction). 

The effect on the signal can be seen in Fig. 5, where the dashed line corresponds to the case when 
the pileup is not included and the continuous line to the one with pileup included. The upper two plots 
show the efficiency of the lepton isolation in the pileup conditions, i.e the fraction of isolated leptons with 
pileup relative to one without pileup as function of the rapidity and transverse momentum of leptons. 
In the lower two plots the effect of pileup on the missing transverse energy and total transverse energy 
is illustrated. Pileup has a small effect on E™ lss distribution as minimum bias events are well balanced. 
On the other hand, it increases the E Ey value by an average value of <~ 800 GeV. 

The most noticeable effect of pileup is the degradation of the lepton isolation. The isolation definition 
includes a rejection at the tracking and at the calorimeter level. For the tracking we require no charged 
particle with a pt > 2 GeV in a cone AR < 0.3 around the electron track (assuming 100 % track-finding 
efficiency). At the calorimeter level we require less than 10 % of the electron energy in cone 0.05 < AR 
< 0.3 around the electron barycenter. On average the isolation efficiency is about 85 % . It is relatively 
independant on the rapidity but increases with p-r to reach a plateau. 
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Figure 5: Pileup effect on isolation efficiency of electrons, on E™ lss and S (solid line corresponds to 
no-pileup case, dashed line - to the pileup included). 

7 Results for electron and muon type leptoquark 

The following results are obtained for muon and electron scalar leptoquarks, Tab. 3. A global 95% 
efficiency for the lepton reconstruction ( ef. ec = 0.95 ) is assumed, based on the results of a detailed 
GEANT study. We deduce a limit of f 470 GeV for a branching ratio of f and a limit of 1200 GeV for 
a branching ratio of 0.5, taking more stringent criterion of the 5 a excess in terms of ^g +B - Fot muon 
type we obtain the same results within less than 10 GeV. In Fig. 6 one can see the lepton-jet invariant 
mass distributions for the Standard Model (SM) background combined with the leptoquark signal for 
various leptoquark masses : 600 GeV, 900 GeV, 1.2 TeV and 1.5 TeV after applying the first level cuts 
and requiring m;; > 150 GeV. The contributions (%) of the various SM background processes after all 
cuts as function of the leptoquark mass are given in Tab. 4, where (a) corresponds to the pr >100 GeV, 
(b): 50 GeV < p T <100 GeV and (c): 20 GeV < p T < 50 GeV. 



M LQ (GeV) 


900 


1200 


1400 


1500 


Signal 


2584 


174.4 


49 


24.6 


Background 


240 


45.27 


11.3 


7.5 


s 

,/S+B 


49 


11.8 


6.3 


4.34 


° = Vb 


167 


25.9 


14.6 


9.0 



Table 3: Number of accepted Ij combinations and significance for L = 100 fb 1 for electron type lepto- 
quark. 
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Figure 6: Invariant lepton-jet mass distributions after first level cuts and mu > 150 GeV for several 
signal leptoquark masses (600, 900, 1200 and 1500 GeV) with expected SM background included. 
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Table 4: Background composition in percents for electron type leptoquark after all cuts 
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8 Mass peak and calibration 



In order to obtain the width we fit the mass distribution for signal only with a gaussian, as the energy res- 
olution dominates over the intrinsic Breit Wigner. The leptoquark mass can be approximately described 
as: 

M LQ = ^(Pi+Pjet) 2 ^/2EiEj{\ - cosB) 

From this relation we can deduce that ^m!^ — "w- " ^ ue to tne excellent lepton energy resolution 
in CMS. Therefore the width of the leptoquark peak is approximately a half of the jet energy resolution. 
For high masses, high energy jets and excluding the Hadron Forward calorimeter (HF) region (\r]\ > 3), 
we expect the constant term to dominate the jet energy resolution. This means that we have (outside of 
cracks and barrel-endcap region): 

zlMml ~ §m = 3.25%. 

MtQ 2 

This is in good agreement with the Fig. 7, where, for example, we obtain a = 42 GeV for Mlq = 1.3 TeV. 

We plot in Fig. 7 the lepton-jet invariant mass applying first level cuts and m;; > 400 GeV, E™ lss < 
180 GeV, £ E T > 1500 GeV, E™ SS /EE T < 0.04. The underestimation of the leptoquark mass is mainly 
due to the finite cone jetfinder and to be corrected by calibration procedure. 
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Figure 7: Leptoquark mass resolution for generated masses of 600, 900, 1200 and 1300 GeV. 



9 Conclusion 

The scalar leptoquark is expected to be observed in CMS detector up to Mlq — 1.5 TeV for 100 fb _1 
and assuming Br(LQ — > e q) = 100 %. The study of the electron and muon leptoquark type gives similar 
results. The width of the leptoquark mass peak depends on the leptoquark mass and is expected to be 
about (7 = 42 GeV at M^q =1.3 TeV. Pileup effect has limited influence, reducing mainly the lepton 
isolation efficiency and shifting SEy value. 
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More details of this study can be found in CMS NOTE available on CMS information server (If 
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